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Abstract 


An approximate polarization correction for ionization loss suggested 
by Sternheivner has been implemented in the galactic cosmic ray trans- 
port code (HZETRN) developed at the Langley Research Center . Sample 
calculations made for the aluminum shield and liquid hydrogen shield 
show no more than a Ai2-percent change in the linear energy transfer 
(LET) distribution for flux compared with those without polarization 
correction. This very small change is expected because the effect of po- 
larization correction on the reduction in stopping power of ions with 
energies above 2 GeV/amu is suppressed by the decrease in galactic 
cosmic ray ion flux at such high energies . 


Introduction 

As the space program proceeds into an era of extended manned space operations, the protection of 
astronauts from galactic cosmic rays (GCR’s) becomes an important problem for mission designers. 
The interaction of high-energy heavy ions originating in deep space with target nuclei results in energy 
degradation and nuclear fragmentations. These fragmentations produce secondary and subsequent- 
generation reaction products that alter the elemental and isotopic composition of the transported 
radiation fields. The transported radiation level and linear energy transfer (LET) distribution inside 
the spacecraft or any critical organ of interest can be estimated by an engineering-design-oriented, 
computer-efficient, galactic cosmic ray transport code (HZETRN) developed at the Langley Research 
Center (ref. 1). Although HZETRN has been widely used, it is continuously being modified for better 
numerical accuracy (ref. 2) and more accurate physics inputs (ref. 3) so that the weight penalty for 
shields can be minimized. 

The stopping power of target materials for high-energy heavy ions is an important input 
physics parameter to HZETRN, besides nuclear fragmentation parameters. At energies above a few 
MeV/amu, electronic stopping power is adequately predicted by Bethe’s theory provided appropriate 
corrections to Bragg’s rule, shell corrections, and an effective charge are included (ref. 4). As the 
input to HZETRN, the high-energy cutoff for the incident GCR spectrum is usually taken to be far 
beyond the minimum stopping power (~2 GeV/amu) where Bethe’s theory starts to overestimate 
and worsens as the energy increases. This overcstimation can be corrected by considering that the 
field of the incoming ion projectile is perturbed by the polarized neighboring atoms. Because the 
incident flux of GCR ions falls off at high energies, Bethe’s theory has been used without correction in 
HZETRN and assumed to be adequate. In this work, an approximate correction for the polarization 
effect (also known as density effect) is included and the effect of this correction is assessed for an 
aluminum shield and a liquid hydrogen shield. 

Galactic Cosmic Ray Transport Method 

When moving through extended matter, heavy ions lose energy through interaction with atomic 
electrons along their trajectories. On occasion, they interact violently with nuclei of the matter 
and produce ion fragments moving in the forward direction and low-energy fragments of the struck 
target nucleus. The transport equations for the short-range target fragments can be solved in closed 
form in terms of collision density (ref. 5). This step leads to the projectile fragment transport as 
the remaining problem of interest. In previous work, the projectile ion fragments were treated as 
if all went straight ahead because of the large kinetic energy possessed by the projectile and the 
fragmentation mechanism. The straight-ahead approximation is found to be quite accurate foi the 
nearly isotropic cosmic ray fluence (ref. 5). 



With the straight-ahead approximation and the target secondary fragments neglected (refs. 5 
and 6), the transport equation may be written as 
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where <j>j{x y E ) is the flux of ions of type j with atomic mass A } and charge Z 3 at x moving along 
the x-axis at energy E in units of McV/amu, a j is the corresponding macroscopic nuclear absorption 
cross section, S(E ) is the stopping power of the protons ( S — S p ), fjk{E, E f ) * s a differential energy 
cross section for the production of ion j in collision by ion /e, and vj is the range scaling parameter 
that is defined as 

v i = j: (2) 

The stopping power S is the sum of nuclear stopping power and electronic stopping power 5 e , 
where S t dominates at high energies. Utilizing the definitions 
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allows equation (1) to be written as 
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which may be rewritten as (ref. 7) 
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where the exponential is the integrating factor with 

C j(r,t) = f ajir + i/jt’) dt' 

J 0 


(7) 


(8) 


Rather simple numerical procedures follow from equation (7). By noting that the first-order 
nature of equation (1) allows , ipj(x 7 r) to be taken as a boundary condition for propagation to larger 
values of x, one may approximate equation (7) as 

rh roc 
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If h is sufficiently small that 
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then, according to perturbation theory (ref. 7), 


tp k {x + h - z,r') « exp[— Ot(r', h - z)] r' + u k (h - z)] 

which may be used to reduce the integral of equation (9), thus resulting in 
ipj{x + h, r) w exp[-< j(r,h)] ipj{x,r + Ujh) 
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Currently, for Zj > 1 and k > j we assume that 

fjk{r, r') = cTjk(r') &{r - r') 

By using equation (13), equation (12) now becomes 
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with r' = r + vjz. Equation (14) is further approximated as 
ipj(x + h, r) « exp[-Cj(r, h)] ipj{x, r + Vjh) 
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Equation (15) is the stepping formalism with energy-dependent cross sections for k > 4 He. The 
corresponding stepping formalism for nucleons has been discussed in detail in reference 4. These 
stepping formalisms arc then used to march the solution from the surface boundary to the desired 
shield thickness. 


Approximate Polarization Correction 

The reduction in the ionization loss of ultrarelativistic charged particles due to the polarization 
of the medium was first treated by Fermi (ref. 8) who assumed that the dispersive properties can 
be described by a single type of dispersion oscillator. Halpern and Hall (ref. 9), Wick (ref. 10), 
and Sternheimer (ref. 11) extended Fermi’s equations to the general case of an arbitrary number of 
dispersion oscillators. In reference 11, Sternheimer evaluated the density-effect (polarization-effect) 
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correction for various elements and compounds and fitted his calculated results by the expressions 


<5 = ° (x<x a ) (16a) 

6 = \nT ] 2 + C + a(xi - x) m {x 0 < x < x x ) (16b) 

6 = hiT ) 2 + C ( x>x x ) (16c) 


where 6 is the density-effect correction that enters the stopping-power formula, 77 = p/m 0 c, p is the 
momentum, m 0 is the rest mass of the charged particle, c is the velocity of light, and 


x = (log 10 e) In 77 = 0.43429 In 77 


The quantities a, m, x 0 , and x x are constants that must be evaluated for each material; C is given 
by 


C — -2 In (I/hi/p) - 1 

where I is the mean ionization potential, hi> p (the plasma energy of the material) is given by 


(17) 


hvp = h{ne* / nm e ) 1 / 2 


(18) 


n is the electron density (in electrons per cubic centimeter), m e is the rest mass, and e is the charge of 
an electron. For compounds or mixtures, the effective mean ionization potential can be determined 

by 
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where n i is the electron density for the ith element and I I is the corresponding atomic ionization 
potential. 


Sternheimer (ref. 11) has suggested that for some practical applications, using only the asymptotic 
density-effect correction (eq. (16c)) may be adequate for all charged-particle energies (ref. 12). 
Armstrong and Alsmiller (ref. 12) subsequently compared the differences in stopping power with 
correct expressions and an asymptotic formula for several elements and compounds. The results 
indicate an overestimate of, at most, 6 percent when using the asymptotic formula. Therefore, only 
the asymptotic formula is used for this study. Bethe’s high-energy approximation to the energy loss 
per unit path, with the density-effect correction 6 included, is given as 
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where N is the number of target molecules per unit volume, Z p is the projectile charge, Z t is the 
number of electrons per target molecule, m is the electron mass, v is the projectile velocity, (3 = v/c , 
and C s is the velocity-dependent shell correction term (ref. 13). Obviously, the use of equation (20) in 
conjunction with the asymptotic formula (eq. (16c)) does not involve any evaluation of the constants 
a, m, and x\ for the material considered. 


Results 

Sample calculations are made for an aluminum shield and a liquid hydrogen shield to explore the 
effect of polarization correction in ionization loss. Figure 1 shows the proton stopping power in these 
two materials predicted by Bethe’s theory with and without Sternheimer’s approximate polarization 
correction. The reduction in stopping power due to the polarization effect occurs at energies beyond 
the minimum value of the stopping power and increases with the energy. However, the effect of 
a large reduction in stopping power on the GCR transport calculation may be suppressed by the 


4 


WFiiffi 




Figure 1. Proton stopping power in aluminum and liquid hydrogen predicted by Bethe’s theory with and without 
Stcrnheimer’s approximate polarization correction. 




5p e , MeV/(g/cm 2 ), in AI 

(c) Iron flux through shield thickness of 15 g/cm 2 . 



(d) Iron flux (£fc through shield thickness of 30 g/cm 2 . 


Figure 2. Comparison of transmitted GCR ion flux through aluminum shield at solar minimum with and without 
Sternheimer’s approximate polarization correction. Flux distribution shown as function of ion stopping power in 
shield. 
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decreasing GCR flux distribution at high energies. Figures 2 5 display the calculated results for 
GCR transport through the aluminum shield (figs. 2 and 4) and the liquid hydrogen shield (figs. 3 
and 5) when using the energy-dependent version of HZETRN. The incident spectrum is based on 
the CREME (cosmic ray effects on microelectronics) model (ref. 14) at solar minimum. For the 
transported proton and iron ions, the flux distributions as a function of ion stopping power arc 
shown with and without the polarization correction in an aluminum medium (fig. 2) and liquid 
hydrogen medium (fig. 3) at a depth of 15 and 30 g/cm 2 . The calculated results based on the use of 
the correction are essentially the same, but, as shown in figure 1, the results begin to diverge beyond 
the minimum stopping power. 
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(b) Proton flux <f> p through shield thickness of 30 g/cm 2 . 



Figure 3. Comparison of transmitted GCR ion flux through liquid hydrogen shield at solar minimum with and 
without Sternheimer’s approximate polarization correction. Flux distribution shown as function of ion stopping 
power in shield. 

Ratios are taken for the LET distributions of the total transported flux through the aluminum 
shield (fig. 4(a)) and liquid hydrogen shield (figs. 5(a) and 5(c)) of various thicknesses to numerically 
quantify such a small difference in transport calculation due to the polarization effect. Deviations 
of, at most, ±2 percent are noticeable even for thicknesses up to 100 g/cm 2 of liquid hydrogen. Most 
of the structured peaks come from the differences in the very high energy region (or very low LET 
region) of each ion species. As more heavy ions are broken up into low LET fragments in penetrating 
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Figure 4. Ratio of differential LET flux calculated with and without Sternheimer’s 
correction. Transmitted GCR flux through aluminum shield thicknesses of 15 and 30 


approximate polarization 
g/cm 2 at solar minimum. 



(a) Total ion flux * through shield thicknesses of 15 (b) Proton flux <f> p through shield thicknesses of 15 and 

and 30 g/cm 2 . 30 s/ cm ' 



(c) Total ion flux <j> through shield thicknesses of 60 (d) Proton flux through shield thicknesses of 60 and 

and 100 g/cm 2 . 100 g/cm 2 . 

Figure 5. Ratio of differential LET flux calculated with and without Sternheimer’s approximate polarization 
correction. Transmitted GCR flux through liquid hydrogen shield at various thicknesses at solar minimum. 


farther into the shield, the low LET ions tend to dominate more strongly. This result is reflected 
in the similarity of structures for the ratios of total flux and proton flux (fig. 4(b) for aluminum 
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and figs. 5(b) and 5(d) for liquid hydrogen) at the low LET region. The humps seen for the liquid 
hydrogen shield at high LET for thicknesses near and beyond 30 g/cm 2 are insignificant because 
most of the heavy ions have disappeared as a result of the large stopping power possessed by liquid 
hydrogen. 

This study concludes that when the full GCR spectrum is considered, the calculated LET spectra 
behind shields are affected by the polarization correction by no more than 2 percent. However, 
the correction is expected to have a larger effect on calculated transport results for the cases with 
geomagnetic cutoff (as in lower Earth orbit) and a much larger effect on laboratory beam experiments 
with energies greater than 2 GeV/amu. 

Concluding Remarks 

An approximate polarization correction for ionization loss suggested by Sternheimer has been 
implemented in the galactic cosmic ray transport code (HZETRN). The effect of a large reduction 
in the stopping powder due to polarization for very high energy ions (beyond 2 GeV/amu) on the 
galactic cosmic ray (GCR) transport calculation is not more than ±2 percent as has been suspected. 
Further detailed implementation with a more accurate formula providing further correction in the 
range from 1 to 10 GeV/amu would be of interest but would provide only a small correction for 
space radiation application beyond the Earth’s magnetosphere. 


NASA Langley Research Center 

Hampton, VA 23681-0001 

February 4, 1993 

References 

1. Wilson, John W.; Chun, Sang Y.; Badavi, Forooz F.; Townsend, Lawrence W.; and Lamkin, Stanley L.: HZETRN: 
A Heavy Ton/Nucleon Transport Code for Space Radiations . NASA TP-3146, 1991. 

2. Shinn, Judy L., and Wilson, John W. : An Efficient HZETRN (A Galactic Cosmic Ray Ti'ansport Code). NASA 
TP-3147, 1992. 

3. Wilson, John W.; Townsend, Lawrence W.; Schimmerling, Walter; Khandelwal, Govind S.; Khan, Ferdous; Nealy, 
John E.; Cucinotta, Francis A.; Simonsen, Lisa C.; Shinn, Judy L.; and Norbury, John W.: Transport Methods 
and Interactions for Space Radiations. NASA RP-1257, 1991. 

4. Wilson, John W.; Townsend, Lawrence W.; Nealy, John E.; Chun, Sang Y.; Hong, B. S.; Buck, Warren W.; 
Lamkin, S. L. } Ganapol, Barry D.; Khan, Ferdous; and Cucinotta, Francis A.: DRYNTRN: A Daryon Transport 
Model. NASA TP-2887, 1989. 

5. Wilson, John W.: Analysis of the Theory of High-Energy Ion Transport. NASA TN D-8381, 1977. 

6. Wilson, John W.: Heavy Ion Transport in the Straight Ahead Approximation. NASA TP-2178, 1983. 

7. Wilson, John W .; and Lamkin, Stanley L.: Perturbation Theory for Charged-Particle Transport in One 
Dimension. NucL ScL & Eng., vol. 57, no. 4, Aug. 1975, pp. 292 299. 

8. Fermi, Enrico: The Ionization Loss of Energy in Gases and in Condensed Materials. Phys. Review, second ser., 
vol. 57, no. 6, Mar. 15, 1940, pp. 485 493. 

9. Halpcrn, Otto; and Hall, Harvey: The Ionization Loss of Energy of Fast Charged Particles in Gases and Condensed 
Bodies. Phys. Review, vol. 73, no. 5, Mar. 1, 1948, pp. 477-486. 

10. Wick, G. C.: On the Retardation of Fast Particles. Nuovo Cimento, ser. 9, vol. 1, June 1943, pp. 302 313. 

11. Sternheimer, R. M.: The Density Effect for the Ionization Loss in Various Materials. Phys. Review , vol. 88, no. 4, 
Nov. 15, 1952, pp. 851 859. 

12. Armstrong, T. W.; and Alsiniller, R. G., Jr.: An Approximate Density-Effect Correction for the Ionization Loss 
of Charged Particles. NucL Instrum. & Methods , vol. 82, 1970, pp. 289-290. 

13. Walske, M. C.; and Bethe, H. A.: Asymptotic Formula for Stopping Power of K-Electrons. Phys. Review , vol. 83, 
1951, pp. 457-458. 

14. Adams, James H., Jr.: Cosmic Ray Effects on Microelectronics, Part IV. NRL Memo. Rep. 5901 (Revised), Naval 
Research Lab., Dec. 31, 1987. 


8 




REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No 0704-0188 

Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, 5uite 1204, Arlington, VA 22202*4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188). Washington, DC 20503 

1. AGENCY USE ONLY (Leave blank) I 2. REPORT DATE I 3. REPORT TYPE AND DATES COVERED 


March 1993 I Technical Memorandum 


4. TITLE AND SUBTITLE 

Polarization Correction for Ionization Loss in a Galactic Cosmic 
Ray Transport Code (HZETRN) 

5. FUNDING NUMBERS 

WU 593-41-21-01 

6. AUTHOR(S) 

Judy L. Shinn, Hamidullah Farhat, Francis F. Badavi, 
and John W. Wilson 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Langley Research Center 
Hampton, VA 23681-0001 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

L-17192 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

NASA TM-4443 

11. SUPPLEMENTARY NOTES 

Shinn and W r ilson: Langley Research Center, Hampton, VA; Farhat: Old D 
Badavi: Christopher Newport University, Newport News, VA. 

dominion University, Norfolk, VA; 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 1 

12b. DISTRIBUTION CODE 


Unclassified Unlimited 


Subject Category 93 


13. ABSTRACT (Maximum 200 words) 

An approximate polarization correction for ionization loss suggested by Sternheimer has been implemented 
in the galactic cosmic ray transport code (IIZFTRN) developed at the Langley Research Center. Sample 
calculations made for the aluminum shield and liquid hydrogen shield show no more than a ±2-percent change 
in the linear energy transfer (LET) distribution for flux compared with those without polarization correction. 
This very small change is expected because the effect of polarization correction on the reduction in stopping 
power of ions with energies above 2 GeV/amu is suppressed by the decrease in galactic cosmic ray ion flux at 
such high energies. 


14. SUBJECT TERMS 

Galactic cosmic ray transport; Stopping power; Polarization effect in ionization loss 

15. NUMBER OF PAGES 

9 




16. PRICE CODE 

A02 

17. SECURITY CLASSIFICATION 

18. SECURITY CLASSIFICATION 

19. SECURITY CLASSIFICATION 

20. LIMITATION 

OF REPORT 

OF THIS PAGE 

OF ABSTRACT 

OF ABSTRACT 

Unclassified 

Unclassified 




1-280 5500 Standard Form 298(Rev. 2-89 


Prescribed by ANSI Std Z39-18 
298-102 


NA5A Langley, 1993 


















